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Abstract 

Within the framework of VT — symmetric quantum mechanics of bound states (which 
works with parity-pseudo-Hermitian Hamiltonians H = VH^V and real spectra) we 
mimic some effects of the double-well structure of potentials by a pair of 5 func- 
tions with mutually complex conjugate strengths. The model is solvable by the 
standard matching technique and exhibits several interesting features. We observe 
an amazingly robust absence of any VT— symmetry breaking. A quasi- degeneracy 
which occurs in the high-energy domain is interpreted as a manifestation of certain 
"quantum beats" . 

PACS 03.65.Fd 03.65. Ca 03.65.Ge 03.65.Bz 



1 Introduction 



In current textbooks, many phenomena (like, e.g., vibrational spectra) and methods 
(like, e.g., perturbation expansions) of quantum mechanics are best illustrated by 
the one-dimensional Schrodinger equation for bound states in a real and symmetric 
well V{x) = V*(x) = V(-x), 



$(x) = Et(x) . (1) 



In such a setting, the transition of the so called VT— symmetric quantum mechanics 
(say, in its form proposed by Bender and Boettcher [1]) may most easily be illustrated 
by an inclusion of an asymmetric and purely imaginary additional potential in the 
same equation, 

V(x) = V s (x)+iV A (x) , V s (x) = V s *(x) = Vs(-x) , V A (x) = V* A {x) = -V A {-x) . 

(2) 

The exactly solvable examples abound [2, 3]. As their most elementary and trans- 
parent sample we may recollect the harmonic oscillator in D dimensions [4] with 



V(x)= : v ~ + 1 l +B(x-ic) 2 , c>0. (3) 
(x — icy 

Here, the centrifugal barrier is regularized and proportional to t = (D — 3)/2 + m in 
the m — th partial wave. 

The studies of models (2) also often involve much more complicated potentials, 
which are not exactly solvable. An encouragement is provided by their more imme- 
diate applicability as well as by the efficiency of many semiclassical [5], perturbative 
[6] and/or purely numerical techniques [7]. In the pioneering paper, Caliceti et al 
[8] paid attention to the very unusual imaginary cubic VT— symmetric anharmonic- 
ity Within sophisticated perturbation theory, she (together with her co-authors) 
succeeded in showing that in spite of the non-Hermiticity of the Hamiltonian, the 
spectrum of energies may remain real. 

Buslaev and Grecchi [9] returned to the more standard, quartic VT— symmetric 
anharmonic oscillator 

V(x) = + \\ + B (x - i cf + D (x - i c) 4 , c>0 (4) 
(x — i cy 
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with a similar motivation stemming from field theory [10]. They were the first who 
demonstrated that the bound-state energies of a non- Hermit ian, VT — symmetric 
model (4) may coincide with the spectrum of a certain Hermitian double-well prob- 
lem. Their sophisticated and explicit, Fourier- type equivalence transformation be- 
tween the Hermitian and non-Hermitian partners may be now better understood in 
the language of A. Mostafazadeh [11]. 

Let us note that the presence of the centrifugal term in eq. (4) played an impor- 
tant technical role in ref. [9], while a similar term has been absent in all the original 
studies of the imaginary cubic forces [5, 8]. Curiously enough, only a re-introduction 
of this term in the generalized model 

V(x) = + B(x-ic) 2 + iC(x-ic) 3 , c>0 (5) 

(x — i c) 2 

seems to have opened the way towards the rigorous proof of the reality of the spec- 
trum of the imaginary cubic oscillators [12]. In the light of this proof, the special 
regular case of eq. (5) with the vanishing £{£ + 1) = is not exceptional at all. 
On the contrary, the inclusion of a strong repulsion £{£ + 1) ^> 1 in eq. (5) seems 
instructive and productive in having paved the way towards the recent clarification 
of the applicability of the current l/£ technique in the VT— symmetric context in 
ref. [13]. Thus, we shall assume that £ is large in all the potentials of the type (3), 
(4) or (5). At this point, we get quite close to the subject of our present short com- 
munication since the implementation of the l/£ technique proved entirely different 
in the Hermitian and VT— symmetric models, and we intend to re-analyze the latter 
case via certain simplified models. 

2 Double wells 

All three VT— symmetric potentials (3) - (5) offer a very good testing ground for 
the comparison of the Hermitian and VT — symmetric calculations. Firstly, they may 
serve as an elementary illustration of the l/£ technique in the Hermitian case because 
their Hermitian versions emerge simply in the limit c — > 0. Secondly, one immediately 
notices that such a step is not mathematically trivial since their centrifugal barriers 
become strongly singular. Fortunately, this does not play any role at all because 
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simultaneously, we have to "shrink" the full real axis of the coordinates to the mere 
positive "radial" half-axis [4]. A fully rigorous discussion of this point is available in 
ref. [9] and enables us to re-interpret our differential Schrodinger equation (1) with 
the real and constrained coordinates x = XHermMan — x (c=o) > as the standard 
central bound-state problem on half-axis. In such a context, the perturbative l/£ 
recipe (as reviewed, e.g., in paper [14] and/or in many other references which are 
quoted therein) is virtually trivial. One simply imagines that the potential V(x) 
possesses an absolute minimum at some real and positive coordinate R > { for 
example, in eq. (3) we get the unique minimum at R = [£(£ + 1)/B} 1 ^, etc }. In 
the vicinity of this minimum, the shape of the potential V(x) may be approximated 
by an appropriate harmonic oscillator well, V(R + £) « const. + oo 2 ^ 2 + 0(1/ R). 
At this point one discovers that the quality of the latter "zero-order" approximation 
increases quickly with the growth of £. Moreover, the expansion of the binding 
energies and wave functions in the powers oil/ R remains feasible and exhibits very 
good convergence properties in practice [14]. 

After a return to the non-Hermitian eq. (1) with x G (—00,00), an application 
of the same idea requires a much more careful analysis. As long as all our sample 
potentials V(x) remain smooth, analytic and confining for all the shift parameters 
c 7^ 0, our wave functions ip(x) stay analytic and normalizable all over the lower half- 
plane of the complex plane of x [1] . At the same time, the manifest non-Hermiticity 
of our sample Hamiltonians at c > leads to the necessity of a re-interpretation of 
their eigenstates. In the spirit of our remarks [15, 16] (cf. also the later reviews [11, 
17]), one must replace the concepts of Hermiticity by pseudo-Hermiticity while the 
unitarity becomes weakened to pseudo-unitarity. As a consequence, the probabilistic 
interpretation of VT— symmetric quantum mechanics acquires the forms which find 
their immediate guidance in the relativistic quantum mechanics, with one of the most 
popular examples provided by the Feshbach-Villars pseudo-Hermitian re-formulation 
of the Klein-Gordon equation [18]. 

The net consequence of all these introductory remarks is that within the frame- 
work of our above exemplification of the VI — symmetric quantum mechanics we are 
free to work with any complex coordinate R of the minimum of V(x). As long as 
we assumed that £(£ + 1) ^> 1 is large, we discover that there exist many minima R 
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generating, in principle, many alternative large— £ expansions. For example, eq. (3) 
with c > gives the minimum of V(x) whenever x A = R A = £(£ + 1)/B. Even in this 
trivial example (where we know all the final solutions in advance!), we have no clear 
criterion for the choice not only between the above-mentioned real R = R(+) = Q 
(with the large g = \[£(£ + 1)/-B] 1//4 | ^> 1) and its negative partner = — g but 
also between the two new purely complex extremes at R — R(i OW er) = —iO and 
R = R{u PP er) = R\iower) = +i @ (with the same real parameter 9 = g^> 1). Moreover, 
the transition to the potentials (4) and (5) (both with B = for simplicity) gives the 
respective rules R 6 = £{£ + l)/(2£>) and i R 5 = 2£{£ + 1)/(3C) so that, generically, 
one has to deal with several pairs of the complex minima such that 

V(R + £)~tt 2 £ 2 , R = R (±) = ± e -i0, VI 2 = nf ±) =uj 2 ±ir] (6) 

where we only know that \R\ ^> 1. In such a setting, we may choose the shift of the 
axis in such a way that c = 9, without causing any change in the spectrum of course. 

The fully exhaustive discussion of the particular cubic example (5) may be found 
in ref. [13]. For the simpler quartic oscillator (4), the same large— £ construction 
would lead to the vanishing and non-vanishing 9 at the positive and negative cou- 
plings D > (occurring, e.g., in refs. [19]) and D < (chosen, e.g., in refs. [9, 20]), 
respectively. Thus, in contrast to the Hermitian case where the real position of the 
absolute minimum of V(x) is unique in the majority of cases of practical interest, 
the regularization mediated by the shifts leads, in accordance with the pattern 
(6), to the most frequent occurrence of the two symmetric minima in V(x) at once. 
Now we may conclude our introductory considerations by a declaration that any 
VT— symmetric double well (6) with the complex (and complex-conjugate) strengths 
Qf±) has not been found solvable in our preceding paper [13], and that it offered the 
main motivation for our forthcoming considerations, therefore. 

3 The model 

We shall try to simulate the effects of the general double attraction in the following 
schematic square-well-like model 

V( x ) = (-u 2 -irf) 5(x + a) + (-U 2 + irf) 5(x - a) x e (-1, 1) (7) 
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where, in the spirit of refs. [3], the forces are reduced to the mere VT— symmetric 
pair of the delta functions at a distance measured by the variable a G (0, 1). This 
means that we plan to solve eq. (1) with the boundary conditions 

V>(±1) = 0, -^ip(±a + 0) - -^ip(±a - 0) = (-lu 2 ± i rj) ip(±a) (8) 

by the standard matching technique as described in textbooks [21]. Under the usual 
"PT— symmetric normalization convention [16] 

il>( x ) = ^s(x)+i^a(x) , ips(x) = ip*s( x ) = ^s(-x) , Mx) = ip* A (x) = -ip A (-x) , 

(9) 

this enables us to denote E = k 2 in the obvious ansatz for the wave functions on the 
three sub-intervals of the whole domain in question, 

= (a — i (3) sin k(x + 1), 
^c( x ) — 7 cosnx + iS sinnx, xe(—a,a) ■ (10) 
^r( x ) = (a + ifi) sin k(— x + 1), 
Moreover, as long as we have 

{ip' L (x) — n(a — % P) cos k(x + 1), 
4>c(x) = — k 7 sin kx + % k 5 cos kx, x G (—a, a) < '< ' : 

ip' R (x) = —k (a + i (3) cos k(— x + 1) 

the appropriate insertions in eq. (8) lead to the four-by-four matrix set of equations 



;i2) 



Here, ji(u) = cosk(1 — a) — uj 2 k~ x sin/t(l — a) and 1/(77) = r\ k^ 1 sin/t(l — a) are 
mere abbreviations. 

We may conclude that the matching conditions may be satisfied if and only if 
the secular determinant F(k) vanishes in eq. (12). This condition has the following 
form 

,2 



x G (—1, —a) 
x G (—a, a) 
x G (a, 1) 

x G (—1, —a) 
x G (—a, a) 
x G (a, 1) 
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sin «(1 — a) 





— cos na 
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F(k) = sin (2k) + — [cos (2k) — cos(2/ta)] + 



u; 4 + rf 

K 2 



sin(2Ka) sin 2 [2/t(l - a)] = . 

(13) 
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Its main merit is its compact form, not quite expected in the light of the previous 
complicated matching formula (12). By its structure it resembles the textbook solu- 
tion of the simple square well so that, in this sense, it is a mere implicit representation 
of the spectrum itself. In fact, one might even believe that an explicit construction 
of this spectrum might be possible (say, in the form of an infinite power series in 
some of the parameters) but due to the utterly elementary character of eq. (13), one 
might still deduce the majority of the relevant features of the spectrum from this 
implicit definition. In order to demonstrate this expectation in detail, the roots of 
eq. (13) will be studied here via the graphs of the function F(k). In this sense, we 
shall determine and display the various qualitative features spectrum of the energies 
E n = K^ n in a series of pictures, via their presentation in the form of the sequences 
of the nodal zeros K n of the function F(k). 

4 Discussion 

The results of our study are sampled in Figures 1 - 7 at the different sets of the 
couplings. In all of these cases we observed, first of all, that the influence of the 
changes of rj is not too relevant so that we have fixed, in all our pictures, rj = 20 for 
the sake of defmiteness. 

In Figure 1 we simulated the system of the two weakly attractive wells which lie 
far from each other. We see that the choice of a = 0.95 and u = 3/2 still leads to the 
mere very weak perturbation of the very well known spectrum of the u = energies 
in square well. The function F(k) oscillates very regularly and its overall amplitude, 
although not quite constant, exhibits just a small slow variation and remains almost 
independent of our choice of the interval of k (the Figure takes k = y/~E G (0, 15)). 

In Figure 2 we see that for the same distance of the wells, the growth of their 
attractive strength up to uo = 15000 leads to the significant growth of the average 
magnitude of the amplitude of F(k) and to the emergence of an "envelope" (i.e., 
an auxiliary curve F av (n) which connects the maxima or minima) with much slower 
oscillations. In addition we discover an overall asymptotic decrease of F av (n) and 
certain less apparent periodic structures when we try to connect some of the "sub- 
dominant" local maxima/minima of F(k). 
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Figures 3 and 4 show the effect of the move of the two wells towards each other. 
Even a not too significant shortening of their distance (to a = 0.85 and a = 0.65, 
respectively) shows that the period of the envelopes becomes shorter and, in the 
latter case, one can see a clear competition between several parallel envelopes which 
start to overlap each other. 

In order to reveal the latter pattern clearly, we omitted the low lying levels. This 
enabled us to see that the interplay between different periods remains fairly regular. 
Step by step, it introduces a much more complicated pattern to the energy spectrum, 
anyhow. The separate levels behave in a more and more chaotic manner, especially 
when we return to the low-lying part of the spectrum in a way illustrated in Figure 
5 where we kept a = 0.65 and weakened u = 150. 

The remarkable feature of all this pattern is that our numerical experiments never 
revealed a merger of the two levels followed, presumably, by their disappearance 
in the complex plane. In the other models, such a phenomenon may exist and is 
usually interpreted as the so called spontaneous VT— symmetry breaking [22]. Here, 
Figure 6 with the fairly small distance a = 0.35 (and with the original very strong 
uo = 15 • 10 3 ) illustrates the situation where one gets very close to such a possibility 
and where several pairs of the real energies appear to be almost degenerate. 

In our last Figure 7 where u = 150, we finally demonstrate that the "chaotic" 
character of the eigenvalues in our schematic and exactly solvable VT— symmetric 
point-interaction double wells may be weakened and partially removed by a return 
to their weaker strengths. We see in the picture that the characteristic pairwise 
"irregular" quasi-degeneracies persist. This type of an irregularity has a slightly 
different source in an ultimate slowdown of the oscillations of the envelopes, the 
"wobbles" of which become comparable to the nodal distances. 

We may summarize that the double-well-like structure of our complex, non- 
Hermitian example complies with some of our intuitive expectations. Thus, the 
"robust" reality of the spectrum or an emergence of the quasi-degeneracy have been 
re-confirmed. At the same time, a few other observations (say, of the multiple and/or 
super-imposed "beats" in the curve F(k)) wait for a deeper understanding and/or a 
generic explanation in the future. 
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Figure captions 

Figure 1. The sinus-like left-hand side curve F(k) of eq. (13) in the weakly perturbed 
square-well regime, i.e., at a "large" distance a = 0.95 of the "shallow" wells with a 
weak strength u = 1.5 and with a (virtually irrelevant) asymmetry measure r) = 20. 
Figure 2. Graphical determination of the high-lying energies E = k 2 (at nodes of 
F(k)) for a very strong attraction at uj — 15 • • • 10 3 (the other parameters are the 
same as in Figure 1). 

Figure 3. A shortening of the "beats" of Figure 2 after the shortening of the distance 
between wells to a = 0.85. 

Figure 4. The emergence of the overlaps between multiple "beats" at a still shorter 
distance a = 0.65. 

Figure 5. A return to the lowest energy levels in Figure 4 and to the weaker attraction 
with uj = 150. 

Figure 6. The quasi- degeneracy of levels for the "deep" wells with uj — 15 • 10 3 at a 
comparatively small distance a = 0.35. 

Figure 7. Same as Figure 6, with a "wobbling-type" irregularity in the energy nodes 
for "shallow" wells with uj = 150. 
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